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K^ , Abstract 

o 

The mixed quark-gluon condensate is another chiral order parani- 
vQ ' etcr in QCD, which plays an important role in the application of 

CP ■ QCD sum rules. In this letter, we study the properties of quark-gluon 

mixed condensate at finite temperature and quark chemical potential 
in the framework of global color symmetry model. Using an infrared- 



O ' dominant model gluon two-point function, we find that the behavior of 

f^ • quark-gluon mixed condensate at finite temperature and chemical po- 

O . tential is similar to that of the quark condensate, and both of them give 

'^ ' the same information about chiral phase transition. We also find that 

the ratio of these two condensates is insensitive to the temperature and 
p\^ • chemical potential, which supports the conclusion obtained recently by 

j^ i the authors using quenched lattice QCD ( They only studied the na- 

ture of the mixed quark-gluon condensate at finite temperature ). 
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The mixed quark-gluon condensates g{q(T^'^G^i,Q) is one of the vacuum 
expectation values which reflect the nonperturbative structure of QCD vac- 
uum. This condensate plays an important role for the low energy phe- 
nomena of hadrons in the framework of QCD sum rules^ [2], where the 
various condensates in the operator product expansion (OPE) are taken as 
input parameters phenomenologically to reproduce various hadronic prop- 
erties systematically. Contrary to the well-known condensates of the lowest 
dimension, (G^^jy) and (qq), the mixed condensate represents the direct cor- 
relation between quarks and gluons in QCD vacuum. Due to the fact that 
the mixed condensate plays a relevant role in OPE as a next-to-leading chiral 
variant operator, the mixed quark-gluon condensate can be used as another 
chiral order parameter of the second lowest dimension for QCD phase tran- 
sition. 

The previous evaluations of the mixed quark-gluon condensate at zero 
temperature and chemical potential were performed in the framework of 
QCD sum rule method 01 HI, quenched lattice QCD|Slin] and QCD ef- 
fective models [3 |Hj, respectively. The value of this condensate is still not 
conclusive from these studies, and the parameter ttIq = g{qcrfj,i,G'^'^q)/{qq) 
determined at the renormalization point u = IGeV ranges from 0.5GeV 
to 3.7GeV^. Recently, the behaviour of the various low-dimensional conden- 
sates at finite temperature and chemical potential is the subject of intensive 
researches due to their direct relation to the QCD phase transition. In con- 
trast to {qq), the thermal properties of the mixed quark-gluon condensate of 
dimension 5 is still least known. Therefore it is interesting to investigate the 
thermal nature of this mixed condensate for it can act as the chiral order pa- 
rameter which is independent of the quark condensation {qq). Furthermore, 
the evaluation of the mixed quark-gloun condensate at finite temperature T 
and quark chemical potential /i has an impact on the hadron phenomenology 
through QCD sum rule. 

So far, the investigation of thermal behavior of g{qa^i,G'^'^q) has been 
performed by the authors [S] through the quenched lattice QCD calculation 
with Kogut-Susskind fermion(KS) field method, where a critical tempera- 
ture Tc ~ 280MeV was obtained for both the chiral symmetry and decon- 
finement restoration transitions. In Ref.[3], the obtained ratio between the 
mixed quark-gluon condensate and the quark condensate is almost indepen- 
dent of the temperature even in the very vicinity of Tc. Since there is sig- 



nificant overlap between contemporary Dyson-Schwinger-Equations (DSEs) 
studies and the numerical simulation of lattice QCD, the DSEs could provide 
an adjunct to lattice QCD. The truncation that is accurate in the common 
domain such as zero temperature and chemical potential and finite temper- 
ature can be used to extrapolate into the finite chemical potential domain, 
which is presently inaccessible in lattice simulations. Therefore, it is neces- 
sary to elucidate the (T, fi) dependence of the mixed quark-gluon condensate 
through DSE-models. In this letter, we will explore the properties of the 
mixed condensate at finite temperature and chemical potential in the frame 
work of the global color symmetry model(GCM). 

As a truncated DSE-model, GCM is a quite successful four-fermion in- 
teraction field theory which can be directly derived through a truncation of 
QCD^Ol- This truncated DSE model has been applied extensively at zero 
temperature and chemical potential to the phenomenology of QCD jlll I12j , 
including the studies of observables from strong interaction to weak inter- 
action area. Furthermore, the truncated DSE models also have made im- 
portant progress in the studies of strong QCD at finite temperature and 
chemical Dotential|13j. 

In a Euclidean space formulation, with {"f^j^u} = '^^nu and 7+ =7^, 
the inverse of dressed-quark propagator at finite temperature and chemical 
potential in the framework of GCM takes the form 

= i^-pA{pk) + i-f4{uJk + il^)C{pk) + B{pk), (2) 

where pk = {p, uj^), ^k = ^k + if^^ ^k = (2^ + l)vrT is the Fermi Matasubara 
frequency, /x is the chemical potential and T,(pk) is the self-energy term of 
the dressed quark propagator. The gap equation for the dressed quark in 
the chiral limit is determined by the rainbow DSE 

{A{pk) - l)f = I / g'D{pk - Qi) ,2r\-^^ri!'^^^2^r,2f~^ (3) 

3 Ji,q A\qi)q^ + C\qi)ujf + B\qi) 

ir<i~ \ TN~2 8 /■ 2 C{qi)pk ■ qi , . 

{C{pk) - 1 Vfc =0 9 D{pk - qi) ^ ^^ ^ 2^ ^ (4) 

1 fi r R /^ ~ ^ 

B{pk) = ^ g^D{pk - qi) ,2/-^^ I r^2^~\-2 , p2r-^ ' (^) 
3 Ji,q A\qi)q^ + C\qi)ujf + B\qi) 

where /^^^ = (3~^ Yl'iL^oo I JTi^ and (3~^ = T. The term g^D{pk) is the 
effective dressed-gluon propagator at finite temperature and chemical po- 



tential, which is treated as an input parameter in the quark gap equation. 
The complex scalar functions A{pi.), B{'Pk) and C{pk) satisfy: F[pk)* = 
-F(p_(fc+i)), F = A,B,C, although not explicitly indicated they are func- 
tions only of (f and tD^. 

We will use the technique proposed in Ref.0 to calculate the mixed quark- 
gluon condensate at finite temperature and chemical potential in the frame- 
work of global color symmetry model. In the context of nonzero tempera- 
ture, the space-time integration f d'^x is replaced by the form f^ dx^ f d^x 
and the continuum integration f d^p in the momentum space is replaced 
by Matasubara frequency summation (3~^ Z]fcL-oo / ^"^P- ^^ ^^^ mean field 
level, the vacuum expectation of any quark operator of the form 

On = (9,,Ag^g,)(g,,Ag^g,) • • . (g,„Ailg.J (6) 

can be calculated straightforward. Here A'*^ is an operator in Dirac, flavor 
and color space. Since the generating functional for the quark field is a typi- 
cal Gaussian type integration at the mean field level in the framework of this 
four-fermion interaction model, the expression for the vacuum expectation 
of quark operator 0„ can be easily derived pi^ 17] and takes the form 

(: 0„ :) = (-1)" E(-)'[4!n ' ' ' 41(^)n^.(i) " " " ('5)...„„,], (7) 

p 

where 11 stands for a permutation of the n indices and S is the dressed- 
quark propagator in the mean field level. For example, we can get the 
familiar expression for the quark condensate {qq)^ at temperature T and 
chemical potential /i 

{qq)t = {-)TrDcS{x,x) 

The trace is to be taken in Dirac and color space, whereas the flavor trace 
has been separated out. The four quark condensate can also be obtained 
easily from Eq.Q and proved to be proportional to the square of the quark 
condensate which is consistent with the vacuum saturation assumption of 
Ref.I?]. 

Due to the glounic degree in giqcr^^G'^'^q), we can not get the expression 
for the mixed quark-gluon condensation directly from Eq. . Since the gen- 
erating functional over gluon field A in this four-fermion interaction model is 



quadratic with a given quark-quark interaction D, the integration over any 
number of gluon fields can be performed analytically [7]. Using the same 
shorthand notation for the typical Gaussian integrations as in Ref. [7j, we 
have 

/"p^e-^^^"'^+^'^ = e-^^""^ (9) 

f VAAe-^^'^^~'^+^^ = {jD)ey^^ (10) 

/p^^^g-iAD-iA+jA ^ [D + {jDf]ey^^, (11) 

where j = q^^—q is the quark color current. It should be noted that D 
stands for the connected gluon two-point Green function at nonzero tem- 
perature and chemical potential and the integration in the exponent takes 
the form f^: dx^ f d^x according to the imaginary time thermal field the- 
ory. Due to the transition from the gluon vacuum average to the quark 
color current ^7/i-2-<? together with the gluon two-point function D which 
is usually treated as an input parameter in truncated DSE-model method 
to get the quark gap equation, we can use the Eq.© to evaluate the mixed 
quark-gluon condensate at temperature T and chemical potential /i in the 
mean field level . Applying the technique described above, we obtain 

g{qix)G^,ix)a'^''qix)) = (-2i)N,^ J [dj^^^g^Diz - x)] 

xtrD[S{z - x)a^uS{x - z)j^] 



{4i)N, f f g^D{zi - x)g^D{z2 

J Z\ J Z2 

-22)7 

(12) 



+(4i)iVe / / g'D(zi - x)g'D{z2 - x 

' Zl J Z2 

xtrD[S{z2 - x)cF^^S{x - zi)7^5(zi - 22)7 



where x = [x^, x) and /^ stands for J^ dx^ f d^x. 
From the dressed quark gap equation in GCM 

T,{x-y) = -g'^D{x-y)j^S{x-y)j^, (13) 

we can get the form 

-g'^D{x - y)S{x - y) = --trz)[S(x - y)] - ---f^trni^ix - y)-f^]. (14) 

Using this equation, the integration over the gluon two-point function D 
in Ea. (|12() is replaced by the quark self-energy S. This technique strongly 



simplifies the calculation of the mixed condensate. The formula for the cal- 
culation of the mixed quark-gluon condensate is now expressed only in terms 
of three scalar functions A,B and C which are solutions of the quark gap 
equation. After Fourier transformation and performing the integration over 
space-time , we get the final integration expression for the mixed condensate 
at temperature T and chemical potential /x. 

B[{2 - A)f + {2 - C))Col] 



g{qG^,a^'''q)i^ = 36 / 

Jk 



fc,, A^q^ + C^Col + B^ 
81 f B^ + 2B[A{A-1)^ + C{C-1)ujI] 



Substituting T = and fi = into Eq. (|T^ . We have 

-A 



g{qG,,a^^''q)lfJ, = (-) ^ !^12 J ds s^ B 



S^2 + 52 



This is just the equation obtained in Ref.[7j. It is clear according to Ea. (|15() 
that in the Wigner phase characterized by i? = 0, the mixed quark-gluon 
condensate takes the value and the nonzero value of this condensate 
signals the nontrivial structure of QCD vacuum. This is the reason that 
giqGfj^^a'^'^q) can also be used as the chiral order parameter as (qq) does. 

It should be noted that our study ignores effects from hard gluonic ra- 
diative corrections to these condensates due to its minor importance for our 
study of nonperturbative properties in the low and medium energy region. 
That means the chosen quark-quark interaction model g^D has a finite range 
in momentum space and the momentum integrals in Ea. (|15() are finite. Be- 
cause GCM is a no-renormalizable effective field theory, the scale at which 
a condensate is defined in our approach is a typical hadronic scale, which is 
implicitly determined by the model gluon propagator g^D and the solutions 
of DSE for the dressed quark gap equations. The similar situation is the 
determination of vacuum condensates within the framework of the instanton 
liquid model|2^ where the scale is set by the inverse inatanton size. 

To explore the properties of the mixed quark-gluon condensate at finite 
temperature and chemical potential, we use the infrared-dominant model 

g^DiPk) = ^{27rfPr^''Sko5ip), (17) 

Id 



which is a generaUzation to T 7^ of the model introduced in Ref.fTT 
where r] ~ l.OGGeV is a mass-scale parameter fixed by fitting tt- and p- 
meson masses. As an infrared enhanced model, it does not represent well 
the behavior oi D^y{pk) away from |p|^-|-a;^ « 0. Consequently some model- 
dependent artifacts arise. However and there is significant advantage in 
its simplicity and since the artifacts are easily identified, the model can 
exhibits many of the qualitative features based on the more sophisticated 
Ansai2;e^lE]- This simple confining model has been used successfully to 
explore the thermal properties of QCDfH], 

Substituting Ea. ((T7|) into Ea. (P3|) . we get an algebraic equation which has 
two qualitatively distinct solutions. The Nambu-Goldstone solution, with 



A(pk) = C(pk 



MpD < 



4 



^(1 -|- ^/l -|- -^2-), otherwise. 



(18) 



r,2 



[ 0, otherwise, 

describes a phase of the model in which: 1) chiral symmetry is dynamically 
broken due to the nonzero quark mass function, B{pi^), in the chiral limit; 
and 2) the dressed-quarks are confined due to the absence of a Lehmann 
representation of the dressed quark propagator. The other one is Wigner 
solution with for which 




B{pk) = 0, A{pk) = C{pk) = -il + Jl + ^), (20) 



which describes a phase of the model with neither dynamical chiral symme- 
try breaking(DCSB) nor confinement. 

In Ref.|14j. this model was used to explore the chiral symmetry and 
deconfinement restoration transition between these two phases character- 
ized by qualitatively different , momentum-dependent modification of the 
quark propagator. The relative stability of the confined and deconfined 
phases is measured by the (T, /x)-dependent vacuum pressure difference or 
bag constant pi7] 

B{T,fi) = P[SNG]-P[Sw], (21) 

where Sng stands for the quark propagator with the Nambu-Goldstone 
solution and S\y means the quark propagator with the Wigner solution. 



P[S] is the pressure obtained by the tree-level auxiliary- field effective action 
in the stationary phase approximation pTj . which takes the form 

P[S] = 'L{TrLn[pS-'] - \Tr[J:S]} (22) 

B(T, fj.) > indicates the stability of the confined phase and hence the phase 
boundary is specified by 

BiT,fi)=0. (23) 

In the chiral limit, the deconfinement and chiral symmetry restoration tran- 
sition in this model are coincident. Substituting Eq. (|18ll2r)|) into Eq. (J23|l . we 
get the eauation[Tl] 

^2 , s&ni^ \2 



^r ■ E r*^l«e(M _ Cfe)-) - ... iS^ll ^ (24) 

^2 A^J^ ^2 ^2 



with ojf < ^ + jU^, p^ = ujf — ojf, which determines the critical line in 
the (T, fi) plane. The critical temperature for zero chemical potential in this 
model is T^ = 0.170Gev which is only 12% larger than the value obtained 
in Ref . |15j using a more sophisticated model and the order of transition is 
the same. This illustrates the simple model's ability to provide a reasonable 
guide to the thermodynamics properties of more sophisticated DSE-models. 
For T = 0, the critical chemical potential is fi^ = 0.3GeV and in case of 
nonzero chemical potential the transition is the first order [nill6j . 

Substituting Eqs. (|18|19|) into Eq.® and Eq. ljT^ . we have the expressions 
for quark condensate and mixed quark-gluon condensate at temperature T 
and chemical potential /i from the infrared-dominant model 

- {M) = ——J2 dpp^MyV^ - 4(p2 - /x2 + ujf + 2ifiuji) (25) 

1 hnax pp 

-g{qaGq) = -^T ^ / dpp^Re{^7j^ - 4(p2 _ ^2 + ^2 + 2ii2u;i)). (26) 



Comparing Eq. (|26j) with Eq. (|25|) . we find that the mixed quark-gluon con- 
densate is proportional to the quark condensate and the ratio between these 
two condensates is a constant independent of temperature T and chemi- 
cal potential fi. It is a nontrivial result that the behavior of the mixed 
guark-gluon condensate at finite temperature and chemical potential is sim- 
ilar to the quark condesate, for g{qaGq) reflects the color-octet component 



of quark-antiquark pairs in the QCD vacuum , while {qq) reflects only the 
color-singlet quark-antiquark components. That means the mixed conden- 
sate g{qaGq) represents the direct correlation between color-octet q — q pairs 
and the gluon field strength G^^, while {qq) represents the direct correla- 
tion between quarks and antiquarks. This nontrivial result was also ob- 
tained in Ref.jS] by quenched lattice QCD simulation with the KS-fermion 
method(note that in Ref.[5], the authors only explored the temperature 
dependence of the mixed condensate and the critical temperature for chi- 
ral symmetry restoration is 280MeV ). The ratio ttiq = g{qcrGq)/{qq) is 
l.QOGeV^ in this infrared-dominant model, which suggests the significance 
of the mixed condensate in OPE. 

In R,ef.|2()j. these two condensates take the form 

giqc^Gq) = i f dX^^^{X\aG\X), (28) 

V J X + rw^ 

where |A) stands for the eigenvector of the Dirac operator as iD\X) = A|A), 
and p{X) the spectral density on A. So the ratio ttiiq can be expressed as 
(A|crG|A)A=o in the chiral limit. The insensitive iT,p) dependence of ttt-q 
suggests that the value of (A|crG|A)A=o has a weak dependence on both the 
temperature and chemical potential in the Nambu-Goldstone phase, even in 
the very vicinity of critical point. 

In Fig.l, the thermal effect on —g{qaGq) is plotted against temperature 
T at several choosen chemical potential points. It is clearly indicated from 
Fig.l that —g(qaGq) decreases continuously to zero with increasing T. The 
chemical potential dependence of —g{qaGq) is shown in Fig. 2 and obviously 
the mixed condensate increases with /_i, up to a critical value Hc{T) when it 
drops discontinuously to zero. This behavior is the same as —{qq), which can 
be attributed to the combination /i^ — tof appearing in Eq.((^SJ) and Ea. H26|) . 

We can see from above analysis that though g{qaGq) and {qq) characterize 
different aspects of QCD vacuum, both of them not only give the same 
critical point for the chiral symmetry restoration, but also have the same 
critical behavior for the QCD phase transition. This conclusion confirms 
the result in Ref.[S] and seems to indicate that all order parameters for a 
phase transition have the universal critical behavior near the critical point. 

In summary, we have estimated the (T, //) dependence of g{qaGq) in the 



framework of GCM. Using a simple, confined DSE-model of QCD, we get 
that the value of the mixed quark-gluon condensate is proportional to the 
quark condensate and the ratio between these two condensates is indepen- 
dent of T and /.i. The large ratio rn^ = l.QOGeV suggests that the mixed 
quark-gluon condensate plays an important role in QCD sum rules. As two 
low-dimensional chiral order parameters, the mixed quark-gluon condensate 
and quark condensate give the same critical behavior for chiral symmetry 
restoration. As mentioned above, this infrared-dominant model is poor to 
represent the ultraviolet behavior of the gluon propagator, so there may ex- 
ist some model-dependent artifacts in our study. Because this model can 
provide a reasonable guide to the thermodynamics properties based on more 
sophisticated DSE-models of QCD, we believe that these results are qualita- 
tively coincident with the more complex Ansdtze . The more sophisticated 
investigation is under progress. 
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Fig.l. The ratio R = {g{qaGq)f^)^^^ /{g{qaGq)f^-ll)''^\ as a function of T for 
a range values of fi. The "structure" in those curves, apparent for small T, is 
an artifact of the fact that the infrared-dominant model does not represent 
well the the quark-quark interaction in the ultraviolet domain [Tl|. 
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Fig.2.The ratio R = {giqo'Gq)^.^^' /{g{qaGq)f^^Q)^5\^ as a function of /i 
for a range values of T. 
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